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This  final  re^x^rt  describes  the  achievements  of  the  AFOSR  sponsored 
research  program  on  laser  interaction  with  dense  metallic  plasmas  under 
Grant  AFOSR-76-2904  during  the  period  1  August,  1975  to  1  October  1979. 
The  principal  aim  of  the  research  effort  is  to  study  the  dynamics  of  a 
high  density,  low  temperature  (10  eV)  plasma  when  illuminated  by  a 
high  intensity  laser  beam.  The  investigation  is  carried  out  both  the¬ 
oretically  and  by  means  of  an  advanced  numerical  simulation  code.  The 
following  sections  give  a  description  of  the  various  efforts  and  accom¬ 
plishments.  A  more  detailed  description  of  much  of  the  work  may  be 
found  in  our  three  Annual  Reports  (1976,  1977,  and  1978)  and  the  several 
Progress  Reports  submitted  eeurlier.  The  first  few  sections  cover  the 
theoretical  effort  followed  by  a  description  of  the  numerical  simulation 
code  investigation. 

A.  STIMULATED  BRILLOUIN  SCATTERING  IN  A  HIGH-DENSITY  2-PINCH  AND 

EXPLODED  WIRE  PLASMAS 

An  important  part  of  the  research  program  concerns  stimulated 
Brillouin  scattering  (SBS)  in  magnetized  plasmas.  Most  important 
effects  are  included  in  our  investigation,  such  as  magnetic  fields, 
plasma  temperature,  arbitrary  angle  of  propagation  of  the  stimulated 
waves,  boundedness  of  the  plasma,  and  density  inhomogeneities.  The 
dispersion  relation  describing  the  stimulated  Brillouin  scattering 
was  derived  and  analyzed. 


The  inclusion  of  a  background  magnetic  field  has  led  to  several 
new  results.  A  magnetosonic  mode  of  the  plasma  is  found  to  be  driven 
across  field  lines  instead  of  the  usual  ion  acoustic  wave  in  the  mag¬ 
netic  field  free  case.  Because  the  damping  rate  of  the  high  frequency 
Alfve'n  wave  is  lower  than  the  damping  rate  of  the  ion  acoustic  wave 
{Te  -  T^)  the  threshold  of  the  instability  is  correspondingly  lowered. 
The  addition  of  the  magnetic  field  inhomogeneity  causes  a  larger  wave 
vector  mismatch  than  in  the  field  free  case.  This  inhibits  growth  of 
SBS  in  nonuniform  plasmas  because  the  waves  convect  more  quickly  out 
of  the  unstable  region  (see  D.J.  Kenney  and  R.S.B.  Ong,  Technical 
Report,  August  1976) . 

B.  FINITE  AMPLITUDE  WAVES  IN  A  MAGNETISED  PLASMA 

The  interaction  of  an  intense  laser  beam  with  a  magnetized  plasma 
generates  density  and  electric  field  fluctuations  in  the  plasma  via  the 
saturation  of  various  instabilities.  Thus  finite  amplitude  waves  are 
present  in  the  disturbed  plasma.  Therefore  it  is  important  to  study 
the  nonlinear  behavior  of  these  finite  amplitude  waves.  In  this  respect 
we  solved  the  problem  of  the  long  time  behavior  of  a  monochromatic 
finite  amplitude  shear  Alfve'^n  wave.  It  is  shown  that  the  nonlinear 
coupling  between  the  Alfven  pump  wave  and  a  background  sound  wave  in 
the  plasma  gives  rise  to  forced  magnetic  sidebands.  The  Alfven  pump 
wave  as  well  as  the  sidebands  steeper,  up  in  the  long  time  scale.  This 
is  shown  by  applying  a  stability  analysis  to  the  nonlinear  SchrSdinger 
equation  derived  for  these  nonlinear  finite  eunplitude  waves  (see 


J.  lonson  and  R.S.B.  Ong,  Technical  Report,  October  1975). 
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C.  DYNAMICS  OF  THE  IRRADIATED  PLASMA:  PONDEROMOTIVE  FORCE 

AND  SOLITON  GENERATION 

We  consider  a  large  an^>litude  electromagnetic  wave  incident  on  the 

critical  surface  (where  the  electron  plasma  frequency  equals  the  inci- 
f 

dent  wave  frequency)  of  a  dense  plasma;  the  ponderomotive  force  generated 
acts  to  expell  matte^'  from  the  region  of  high  field  intensity  forming  a 
density  depression.  This  depression  serves  to  trap  the  electric  field, 
enhemcing  the  ponderomotive  force,  which  in  turn  deepens  the  cavity. 

In  this  manner  the  trapped  electric  is  able  to  reach  a  very  large 
amplitude.  Its  envelope  has  the  shape  of  a  solitary  wave  cuid  is  des¬ 
cribed  by  a  Schrodinger  equation  with  a  nonlinear  potential.  The  non¬ 
linearity  in  the  potential  is  of  the  exponential  type;  it  shows  an 
exponential  decrease  as  the  electric  field  increases.  In  many  treat¬ 
ments  of  the  problem  the  exponential  term  is  usually  expanded  ^md 
higher  order  terms  in  the  series  expansion  are  then  neglected.  The 
consequence  of  such  a  truncation  is  that  the  density  becomes  unphysical 
as  the  electric  field  grows  instead  of  decreasing.  On  the  other  hand 
the  inclusion  of  the  exponential  forces  the  field  to  saturate  at  large 
values  in  agreement  with  observation. 

We  have  retained  the  exponential  form  of  the  potential  in  the 
nonlinear  SchrSdinger  equation  and  solved  for  the  envelope  numerically. 
The  result  is  a  solitary  pulse  or  "soliton".  We  also  showed  that  the 
soliton  solution  is  a  singular  solution  of  the  nonlinear  Schrodinger 
equation  possessing  an  exponential  nonlinearity. 
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D.  THERMAL  SELF-FOCUSING  OF  LASER  LIGHT  IN  A  DENSE  PLASMA 

The  problem  investigated  in  this  part  of  the  research  arises  from 
experimental  results  observed  by  the  experimental  group  in  our  labora¬ 
tory.  A  CO^  laser  is  directed  perpendicularly  to  an  over-dense  plasma 
column,  and  the  transitted  intensity  of  the  beeun  is  measured  as  a  func¬ 
tion  of  the  incident  power  level.  Two  results  have  been  observed: 
first,  the  transmitted  intensity  shows  a  threshold  effect  as  a  function 
of  the  incident  power.  As  the  incident  power  is  increased,  the  trans¬ 
mitted  intensity  jumps  discontinuously .  Below  a  certain  laser  intensity 
no  transmission  is  observed.  Second,  when  a  photographic  plate  is 
placed  to  observe  the  transmitted  light,  a  Fresnel  diffraction  pattern 
is  recorded.  This  pattern  suggests  that  the  trcinsmitted  beam  emerges 
from  a  sharp  hole  in  the  over-dense  plasma. 

In  order  to  explain  the  above  results  we  make  use  of  the  theory  of 
thermal  self- focusing  of  a  laser  beam  in  a  plasma.  As  the  laser  beam 
interacts  with  the  plasma  the  large  amplitude  oscillating  electric 
field  of  the  beam  couples  to  the  electrons  in  the  plasma.  The  fre¬ 
quency  of  the  oscillating  wave  electric  field  is  too  high  for  the 
massive  ions  to  follow.  The  process  transfers  electric  field  energy  to 
kinetic  energy  in  the  electrons.  The  net  effect  is  an  increase  in  the 
electron  temperature  inside  the  beaun  area.  The  energy  of  the  electrons 
acquired  from  the  incident  laser  beaun  is  calculated  as  a  function  of 
the  laser  intensity.  From  the  plasma  parameters  in  our  laboratory 
experiment  (cf.  Annual  Report  III)  we  determine  that  the  main  energy 
dissipation  mechanism  is  electron  heat  conduction  in  the  direction 
tramsverse  to  the  axis  of  the  incident  laser  beaun. 
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Next  the  electron  density  distribution  in  the  tremsverse  direction 
IS  calculated  from  the  steady  state  energy  balcince  equation.  From  the 
value  of  the  sound  speed  in  our  plasma  laboratory  we  determine 
that  the  ch^u■acteristic  time  scale  for  hydrodynamic  motions  transverse 
to  the  beam  is  larger  than  the  tender ature  equilibration  time.  Thus 
we  may  neglect  hydrodynamic  motions  in  the  tr^ulsverse  direction.  This 
implies  that  the  electron  pressure  is  approximately  constant  across 
the  laser  beam.  As  a  result  the  transverse  electron  density  is  inversely 
proportional  to  the  electron  temperature  distribution.  Thus,  for  an 
incident  laser  beam  with  an  electric  field  having  a  Gaussian  profile  the 
electron  heating  is  nonuniform  in  the  transverse  direction.  This  means 
that  the  electron  temperature  is  maximum  on  the  be2un  axis  and  decreases 
in  the  direction  perpendicular  to  the  axis.  Consequently  the  electron 
density  is  minimum  on  the  axis  of  the  incident  laser  becun.  This 
depression  in  the  electron  density  increases  the  dielectric  const^mt 
(and  also  the  index  of  refraction)  of  the  plasma.  It  is  larger  inside 
the  beam  than  outside.  Since  light  rays  bend  towards  regions  of  higher 
dielectric  constant,  the  beeim  focusses  in  towards  the  centerline  area. 
This  causes  a  further  increase  in  the  field  intensity  in  the  region  near 
the  beam  axis  with  a  corresponding  increase  in  the  value  of  the  dielec¬ 
tric  constant  of  the  plasma.  This  is  the  process  of  thermal  self- 
focusing  and  it  enables  the  laser  beam  to  tunnel  itself  through  ai\ 
over-dense  plasma. 

We  have  mathematically  analyzed  the  thermal  self- focusing  of  the 
laser  beam  in  our  laboratory  experiment.  The  threshold  is  determined 
and  the  laser  amplitude  profile  calculated  as  it  propagates  through 
the  plasma  (see  A.  Schmutt  and  R.S.  B.  Ong,  Technical  Report, 


February  1979) . 
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E.  DYNAMICS  OF  THE  STRONGLY  TURBULENT  PLASMA  GENERATED  BY 

LASER-PLASMA  INTERACTIONS 

A  high  intensity  laser  beam  disturbs  a  plasma  very  violently  such 
that  spatially  localized  electric  fields  (solitons)  may  appear.  These 
localized  intense  fields  are  basic  in  the  understanding  of  strongly 
turbulent  plasmas.  Further,  such  localized  nonlinear  entities  can 
interact  strongly  with  individual  electrons  and  ions  thus  leading  to 
a  signific^mt  modification  of  the  zero-order  properties  of  the  plasma 
itself  (e.g.  density,  temperature,  flow  velocity,  etc.). 

It  is  recognized  that  these  localized  intense  fields  (solitons) 
in  unmagnetized  plasmas  are  subject  to  a  collapsing  type  of  instability 
in  more  than  one  dimension.  Of  course  in  a  real  plasma  various  mecham- 
isms  will  intervene  to  prevent  total  collapse.  In  particular,  if  the 
nonlinear  term  is  the  full  exponential  one,  the  soliton  is  stable 
against  collapse.  However,  this  condition  implies  that  the  plasmon 
energy  density  is  of  the  order  of  the  thermal  energy  density.  But  it 
is  well  )cnown  that  such  high  levels  of  turbulence  are  not  attainable  in 
plasmas  since  various  mechamisms,  such  as  nonlinear  Landau  damiping, 
intervene  to  transfer  energy  from  the  waves  to  the  random  particle 
motions. 

The  presence  of  a  large  bacltground  magnetic  field  can  improve  the 
situation.  It  is  with  this  regard  that  we  investigate  whether  a  laurge 
bac]cground  magnetic  field  could  provide  the  rigidity  which  would  pre¬ 
vent  the  collapse  of  solitary  waves. 

We  consider  finite  amplitude  solitary  Alfven  wave  in  a  collision¬ 
less  plasma  with  a  large,  externally  applied  magnetic  field.  Recent 
studies  showed  that  in  the  absence  of  density  perturbations  (i.e.,  the 
incompressible  MHD  case)  nonlinear  Alfven  waves  form  an  axially  symmetric 


w.tveijuiJe  jIoiw  tJie  ambient  mavjnetic  fielJ.  The  Alfven  wave  exhibits 

a  vortex-liK.e  behavior  insiilo  the  wavevjuide.  These  results  were 

obtainevi  by  retainiiivi  second  order  nonlinear  terms  in  the  HMD  e^iuations 

describing  Alt'ven  waves  in  cylindrical  geometry.  In  addition,  exact 

solitary  Alfven  wave  solutions  of  arbitrar%'  amplitudes  in  two- 

dimensional  pl.ine  'jeometr^'  are  obtained  by  Haseoawa  and  Mima  based  on 

Kinetic  theory.  In  their  worK  corrections  to  all  orders  in  the 

■implitude  consistent  with  the  low  3  =  nT  (B  '  B’T)  approximations 

e  o 

are  taken  into  account.  The  nonlinear  wave  is  shown  to  have  a  loncj 
structure  alonu  the  ambient  field  line.  Furthermore,  transverse 


dispersion  .ind  nonlinearity  uive  rise  to  a  localization  in  the  per¬ 
pendicular  direction  with  a  characteristic  lencjth  of  an  ion  Larmor 
radius.  In  t.his  pap'er  we  present  a  nonlinear  Kinetic  Alfve^n  wave 
structure  with  axial  symmetrs'  about  the  line  of  force  of  the  bacKaround 
maubetic  field.  The  full  nonlinearity  due  to  the  coupling  of  the  den¬ 
sity  with  the  longitudinal  electric  potential  is  taken  into  account 
to  all  orders  in  amplitude  consistent  with  the  low  3  approximation. 

We  assume  a  low  3  plasma  with  a  constant  background  magnetic  field 
in  the  z-direction.  The  low  3  assumption  enables  us  to  use  the  two 
potential  fields  and  V  to  describe  the  electric  field,  •  -  34  '3z 

and  E  «  -  7  ^  ■  -  r3^ '3r  ,  where  we  have  assumed  axial  <\Tnmetr\' 

(  3/33  •  d) .  These  fields  produce  only  shear  perturbations  in  the 


magnetic  field.  Hence  we  have  B 


& 

"■br  32  -J  "  3r3z 


con St. ant ,  B  -  d ,  and 
r 


(d  -  V) 


where  E  is  the  radial  electric  fielvl. 
r 
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Ampere's  law  yields  the  following  relations: 

.  -  «  ■ 

32  c 

-  :r-  (rB„)  =  —  : 

r  3r  o  c  z 

Because  of  the  low  0  assumption  the  dominant  contribution  to  the 

longitudinal  current  is  made  by  the  electrons.  Hence  we  may  write 

3j  3(e  n  ) 

z  e 


(2) 


(3) 


(4) 


3z  3t 

The  low  frequency  characteristic  of  the  Alfven  wave  motion  enaibles 

the  electrons  to  thermalize  along  the  magnetic  field  lines.  Thus  the 

electron  density  couples  to  the  longitudinal  electric  potential  ^  via 

the  Boltzmann  distribution 

n  =  n  exp  (ei^z/T  )  (5) 

e  o  e 

The  following  transformations  are  now  introduced  in  order  to  make 
the  system  of  Eqs.  (l)-(5)  dimensionless: 

P  =  rQ/c 


c  =  /T  /m. 
s  e  1 


z  =  z(a)  ,/c)  =  z(a)/V  ) 
pi  A 


=  e  B  /m .  c  ; 

o  1 

2  *5 

(ij  .  =  { 47:n  e  /m . ) 
pi  o  1 


T  =  Qt 


ij)  =  et|;/T 


=  e4i/T 


E  =  E  ec/O)  .  T 
z  z  pi  e 


E  =  E  e  c  T 
r  r  s  e 


;  n  =  n  /n 
e  o 

B.  =  Ba/(47rn  T 
t)  0  o  e 


where  n^  is  the  unperturbed  plasma  density.  Equation  (1)  becomes 


3B-  3e  3e 

_ 9  ^  _ z  r 

3T  9p  9z 


3^ 


(0  -  ij)) 


(6) 


Combination  of  Eqs.  (3)  and  (4)  yields 

3B 


.-1  3 


3  f  o  Q  I  _  ^ 


3p 


3z 


3n 

3T 


(7) 
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In  the  low  (i  approximation  the  electrons  are  tightly  bound  to  the 
magnetic  field  lines.  Also  the  longitudinal  ion  current  may  be  neg¬ 
lected.  The  ion  current  is  primarily  due  to  the  polarization  drift 
in  the  radial  direction.  This  is  obtained  by  considering  the  ion 
guiding  center  acted  upon  by  an  oscillating  electric  field  perpen¬ 
dicular  to  the  background  magnetic  field.  For  low  frequencies  (U)  <<  Q) 
the  ion  equation  of  motion  yields  an  expression  for  the  ion  polariza¬ 
tion  drift 


3e 


where  is  the  ion  velocity  in  the  radial  direction  and  is  the 
radial  electric  field. 

Substituting  this  into  Eq.  (2)  we  obtain,  in  dimensionless  units, 


SB.  3e 

_ 9.  _  ~  _ r 

3z  3t 


(8) 


Equation  (8)  shows  the  coupling  of  the  B  -field  to  the  transverse 

y 

electric  field  via  the  ion  polarization  drift. 

If  we  choose  a  coordinate  axis  5  -  K  z  -  T  moving  along  z, 

z 

and  scale  DtoK,i.e.  ,p  =  Kp,  where  K  =  k  cU/u)  OJ  .  =  V,/(w/k  ) , 
^z  '^z'^  zz  piAz 

we  may  integrate  Eq.  (7)  and  obtain 


“6> 


(9) 


We  rewrite  Eq.  (6)  as 


3B  3E,  3E^ 

- i  _ _ 5.  -  K  — 

3C  3p  2  3^ 


and  Eq.  (8)  as 


3B. 


2  3C 
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Eliminatinq  (3E^/3?i) ,  we  obtain 


‘*^1  ■  3C  ■  "  3o 


With  the  aid  of  Eq.  (5)  we  may  write 

'  ~-l  3n 

z  z  3i, 

Now,  eliminatinq  between  the  last  two  equations  results,  after  a 
little  manipulation,  in  the  following  relation 

3B„ 


'  ■  ■"z  5F  JK 


3  3n  ~-l  3n  3n 

-  K  _  ^  ^  ^  n  ^  ^ 


We  now  use  Eq.  (9) ,  and  define  the  operator 


^  3^7  i  ^  ^ 


such  that ,  e.g. , 


If  -^)5f 


We  apply  these  notations  to  our  equation  for  (38^/3^)  above  and 
obtain,  after  some  calculations,  the  following  nonlinear  wave 
equation  in  : 


3B 


3B 


3B, 


"'V  -  ST  -  -ilsr  *1"“  |f(‘’  --7)“e 

Let  US  assume  B  to  be  harmonic  in  i.e.  B  »  f (tc^)  exp  (itc^), 

•> 

whore  <“  ■  K^*"  -  1.  Making  use  of  Eq.  (9)  ,  the  wave  equation,  (10)  , 


then  becomes 


(--A) 


B  -  n  (K  “  -  n)  B,  -  0 
0  z  0 


Define  F  =  1  -  n.  Then  using  K  “  ■  1  +  tc  ,  we  can  write  the  last 
equation  as 
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Bq  -  (1  -  F) {K  +  F)  «  0 


or 


M  'k  -  (1  -  <^)  F  -  <J 


Bg  -  0 


(11) 


o(?5-L  ^  ^ 

Away  from  the  ambient  magnetic  field  line,  i.e.,  when  P  is  large,  the 
solution  of  Eq.  (11)  may  be  approximated  by  a  modified  Bessel  function 
K^(KP).  Close  to  the  magnetic  field  line,  P  «  1,  the  near  field  solu¬ 
tion  is  approximated  by  Ij^(icp).  Connecting  these  near  and  feu:  field 


approximations,  in  the  region  where  P  -  1,  is  a  solution  of  the  form 

of  a  Bessel  function  of  the  first  Itind,  Ji(1f„^  -  (1  -  K^)  F  - 

where  F„  is  some  mciximum  value  of  F.  Equation  (11)  has  been  solved 
M 

numerically  (Table  1)  and  Figure  1  shows  the  solution  for  various 
values  of  the  parameter 


S  r 


P) 


(Table  1) 

(Figure  1) 

(Figure  2) 

The  B  -field  is  localized  transverse  to  the  bac)cground  magnetic  field 
b 

line  with  a  characteristic  spatial  dimension  of  the  order  of  a  few  ion 
Larmor  radii.  Thus  it  forms  a  waveguide  along  the  ambient  magnetic 
field  line.  By  inspection  of  the  far  field  solution,  K^(Kp),  we  may 
determine  the  width  of  the  solitary  wave,  p  .  The  value  of  the  argu- 
ment  where  the  amplitude  goes  to  zero,  is  given  by. 


(12) 


from  which  we  obtain 


p  =  3K  /(K  ^  -  I)*’ 
M  Z  Z 


(13) 


Figure  2  shows  how  the  width  of  the  solitary  wave  varies  with  K^. 

For  small  an^litude  solitary  waves,  is  somewhat  larger  than  unity, 

cuid  the  denominator  in  Eg.  (13)  is  small.  Thus  the  solitary  wave 

spreads  (see  Fig.  1).  However,  as  increases,  the  denominator  in 

Eq.  (13)  tends  tow2u:d  K  ,  and  the  wave  contracts.  Thus  K  parametrizes 

z  z 

the  nonlineeurity.  The  width  contracts  to  its  asymptotic  value  of 

approximately  three  ion  Larmor  radii.  The  eunplitude  of  the  Bg-field 

peaJcs  around  p~l,  and  looking  along  the  ambient  magnetic  field  line 

the  solitary  wave  appears  as  a  ring-shaped  region  of  large  magnetic 

field.  This  is  shown  in  Fig.  3  where  the  amplitude  of  the  azimuthal 

magnetic  field  B  is  shown  proportional  to  the  density  of  the  circles. 

2 

In  the  opposite  limit,  as  K  (euid  thus  the  nonlinearity)  goes  to  zero, 
the  soliteUiY  wave  spreads  to  infinite  width  destroying  localization, 
and  yielding  a  plane  wave-type  solution.  As  previously  shown  we  may 
re-write  Eq.  (11)  as 

(v  ~  ^)  -  n)  Bg  =  0 

2 

from  which  it  is  clear  that  as  n  eind  approach  unity,  we  are  left 
with  a  linear  equation. 


(Figure  3) 


1.0,  0.7,  AND  0.5  CALCULATED  ON  THE 


TABLE  I. 


Figure  1 


Figure  2 


Figure  3 


VALUES  OF  f  FOR  </K  = 

z 

UNIVERSITY  OF  MICHIGAN  AMDAHL  470  V/6. 

.  Radial  Profile  of  the  Azimuthal  Amplitude  f,  for  Some 
Values  of  where  Bg  =•  fCicp)  exp  (iK^)  . 

2 

.  Solitary  Wave  Width  as  a  Function  of  K 

M  z 

.  Azimuthal  Magnetic  Field  Amplitude  in  the  p  -  0  Plane 
Perpendicular  to  the  z-axis.  The  Ambient  Field 
B^  =  (0,0, B^)  is  Directed  into  the  Page.  The  Amplitude 
of  the  Wave  Magnetic  Field  B  =  (0,Bg  (icp)  ,  0)  is  Shown 
Proportional  to  the  Density  of  the  Circles. 
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k/K  =  1.0 

k/K  =  1.0 
'  z 

ic/K  =0.7 
z 

k/K  =0.5 
z 

P 

f 

f 

f 

0.00 

0.000 

0.000 

0.000 

0.  25 

0.  108 

0.048 

0.023 

0.  50 

0.  225 

0.  102 

0.049 

0.75 

0.350 

0.  154 

0.082 

1.00 

0.466 

0.210 

0.111 

1.25 

0.499 

0.  258 

0.141 

1.  50 

0.431 

0.312 

0.173 

1.75 

0.298 

0.321 

0.210 

2.00 

0.201 

0.  275 

0.229 

2.25 

0. 137 

0.223 

0.248 

2.  50 

0.088 

0. 181 

0.249 

2.75 

0.042 

0.  141 

0.237 

3.00 

0.002 

0. 108 

0.212 

3.  25 

0.000 

0.077 

0.183 

3.  50 

0.000 

0.054 

0. 152 

3.75 

0.000 

0.037 

0. 126 

4.00 

- 

0.022 

0. 106 

4.25 

- 

0.011 

0.087 

4.  50 

- 

0.000 

0.069 

4.75 

- 

0.000 

0.054 

5.00 

- 

0.000 

0.041 

5.  25 

- 

- 

0.031 

5.  50 

- 

- 

0.019 

5.75 

- 

- 

0.008 

6.00 

0.000 

0.000 

0.000 

Table  1.  Values  of  f  for  k/K^  =  1.0,  0.7,  and  0.  5  calculated  on  the 


U.  of  Mich.  Amdahl  470  V/6. 
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F.  C0U.1SI0NAL  EFFECTS  ON  WAVT-BREAKING  IN  A  COLD  PLASMA 

As  the  laser  beam  in  our  laboratory’  experiment  tunnels  its  way 

through  the  o\’er-dense  plasma  a  large  density  gradient  transverse  to 

the  beam  is  generated.  As  the  plasma  is  owr  dense,  there  exists  a 

region  of  critical  densitv  at  some  radius  r,,  where  0  <  r.  <  r  ,  and  r 

1  1  o  o 

13  the  edge  of  the  incident  laser  beam.  In  the  neighborhood  of  r^^, 
the  region  of  critical  density,  resonant  absorption  of  the  laser  beam 
takes  place.  It  is  well-known  in  collisionless  plasmas  that  resonant 
absorption  is  one  of  the  effective  ways  in  which  an  electromagnetic  wave 
can  deposit  large  amounts  of  energy'.  In  the  process  a  large  amplitude 
electrostatic  wave  is  generated  which  e\’entually  breaks  with  the  conse¬ 
quent  production  of  energetic  isuprathernal)  electrons.  Ion  dynamics 
and  thermal  effects  modify  the  wa-e  breaking  eunplitude.  They  reduce 
the  role  of  the  wave-breaking  process  of  energy  dissipation  by  provid¬ 
ing  competing  modes  of  energy  release.  Collisions  produce  a  similar 
effect.  As  our  laboratory  plasma  is  i-ery  dense  and  has  a  relatively 
low  temperature,  collisional  effects  are  important.  Therefore  u’e 
investigate  its  effect  on  the  wave-breaking  amplitude. 

Consider  an  electromagnetic  wave  incident  upon  a  plasma  with  a 
linear  density  profile.  Let  E  be  the  component  of  the  incident  wave 

di 

electric  field  parallel  to  the  density  gradient.  A  Lagrangian  formalism 

will  be  used  to  describe  the  nonlinear  electron  plasma  oscillations  under 

the  influence  of  the  driver  field  E  .  From  the  equation  of  motion  of  the 

d 

electrons  (including  the  collision  term)  and  Poisson's  equation  it  is 
found  that  the  electron  displacement  >S  satisfies 

Vi5  "  vx  )  5  -  t'  sin  t  (F-l> 

pe  o 


where  t 


-  e  E,/m  L;  u),  is  the  frequency  of  the  incident  electro- 

d  d  u  d 

maqnetic  wave;  L  is  a  character  istic  distemce,  e.g. ,  the  density 
gradient  scale  length.  All  distances  and  times  in  Eq.  (F-1)  are 
normalized  with  respect  to  L  and  u)^  ^  respectively.  Thus  Eq.  (F-1) 
is  an  equation  in  dimensionless  form.  The  effect  of  collisions  is 
taken  into  account  by  the  second  term  in  Eq.  (F-1)  where  v  is  propor¬ 
tional  to  the  electron- ion  collision  frequency.  The  displacement  of 

an  electron  initially  at  position  x  is  given  by  S,  i.e. 

o 

iS(x  .t)  -  x(t)  -  X  ,  (F-:) 

o  o 

It  has  been  implicitly  assumed  in  Eq.  (F-1)  that  5  <<  1. 

The  solution  of  Eq.  (F-1)  is  given  by  (see  Appendix) 


6(t)  »  [e/cos  Y  (uip"  -  1  +  V  tan  Y)  1  (exp  (-  vt/I)  ] 

.  ‘j. 


j(sin  y)cos  [(a)  ‘  -  i  v‘  )tj 


V  sin  Y  -  cos  y)  sinf  /  u)  ^  -  -7  v"}  1 1  //u)  “■  -  7  v‘  )  *  ) 

2  ^'pq  J 


sin  (t  -  y) 


(F-3) 


where  Y  =  tan  (v/uip"  -  1).  Let  us  now  ccnsidar  the  behavior  of  ^(t) 
at  the  resonance  region  where  the  electrostatic  wave  grows  most 


rapidly.  We  take  the  limit  as  1  and  let  v*"  <<  uj  *.  The  result 

P  P 


is 


5(t)  -  (e/v)  exp  (-  vt/2)  ^  cos  t  +  7V  sin  t )  -  cos  t^  (F-4) 


Note  the  limiting  effect  that  the  exp  (-  vt,  2)  factor  has  on  the  wave 
amplitude.  Clearly,  for  sufficiently  large  values  of  v  the  electro¬ 
static  wave  is  drastically  modified.  It  is  also  to  be  obser\*ed  that 
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in  the  limit  V  •*  0  the  collisionless  result  of  Koch  and  Albritton  is 
recovered. 

In  order  to  examine  the  effect  of  the  collision  frequency  in  more 
detail  we  consider  the  time  t^  at  which  the  wave  breaks.  This  condi¬ 
tion  occurs  when  3i5/3x^  -  -  1,  i.e.,  when  the  electron  paths  cross  one 

anot.her.  We  now  apply  this  condition  to  Eq.  (F--J)  and  find  that  t, 

b 

satisfies  the  relation 

2  1 
v'/e  »  1  -  exp  (-  vt^/2)(l  +  —  vt^)  .  (F-5) 

This  relation  is  shown  graphically  in  Fig.  1.  The  important  thing  to 
note  is  that  the  right  hand  side  is  never  larger  than  unity.  Thus, 
if  V  >  the  electrostatic  wave  cannot  break.  Instead  of  the  wav’e 

to  break  and  disperses  its  energy  i.nto  hot  electrons  that  stream  out 
of  the  region,  the  electron-ion  collisions  thermalize  the  wave  energy 
in  the  resonance  region. 

Thus  to  prevent  wave-breaking  the  electron-ion  collision  frequency 

2  k 

has  to  exceed  (e  E,/m  L)  .  For  lower  values  of  v  the  wave  does 
d  e  d 

break,  but  is  now  modified  by  collisions.  Aside  from  the  much  longer 

breaking  time,  the  resulting  electron  velocity  is  also  reduced. 

The  electron  velocity  at  wave-breaking  V  is  obtained  by  taking  the 

c 

time-derivative  of  5  in  (F-4)  and  evaluating  it  at  t  =  t  .  Figure  2 

b 

shows  the  electron  velocity  at  wave-breaking  as  a  function  of  the 
collision  frequency.  Thus  collisions  can  significantly  reduce  the 
energy  carried  away  by  the  hot  electrons. 

Finally,  the  limitations  of  this  analysis  need  to  be  pointed  out. 
With  the  longer  wave-breaking  times  the  effect  of  ion  bunching  becomes 
more  important.  However,  this  can  be  offset  somewhat  by  consideration 
of  ion  collisions.  Furthermore,  as  collisions  have  the  effect  of 


themaliiincj  the  wave 


enervjy  in  the  resonance  reijion, 
would  also  become  significant.  Finally,  it  is  tacitly 
above  analysis  that  vt  -  1,  i.e.,  the  electrons  suffer 
collision  before  the  wave  breaks. 


temperature  effects 
implied  in  the 
at  least  one 


(Fi^gure  ) 


(Figure  ) 


APPENDIX  (to  Section  F) 


Consider  the  equation 


6+v6+a)  (x)5=e  sin  t 

pe  o 


(F-A.l) 


The  general  solution  of  the  associated  homogeneous  equation  is, 
2^1  2 

assuming  ^  ■4  ^  • 

Oj,  =  exp  I  -  7Vt)[c^  +  c^  sini 

with  Cj^  and  c^  to  be  determined  from  the  initial  conditions. 

To  find  a  particular  solution  of  the  nonhomogeneous  equation, 
(F-A.l) ,  we  try 

6p  =  c^  sin  (t  -  Y) 

Substituting  this  into  (F-A.l)  and  setting  the  coefficients  of  the 
sin  t  emd  cos  t  terms  equal  to  zero  we  obtain 

Y  =  tan  ^  [v/(w  ^  -  Ij] 

c 

2 

c^  =  e/cos  Y  (<*Jp  -  1  +  V  tcui  Y) 

Thus  the  solution  so  far  is 

6(t)  =  (c^  cos  at  +  c^  sin  at)  exp  |  j  vtj  +  c^  sin  (t  -  y) 


2  X  2  ~ 

with  c-  and  y  defined  above  and  a  =  foj  -  -7  .  We  now  determine 

3  '•  p  4  ^ 


c^  and  c^  from  the  initial  conditions 
5(t=0)  =  0  , 


6(t=0)  =  0 


The  final  form  of  the  solution  is  then  given  by 


l^/cos  Y  -  1  +  V  tcin  Y)j  ^ 

,  I  exp  (-  vt/2)  [sin  y  cos  |  ~  ^  t  ^ 


+  ( j  V  sin  Y  -  cos  Y  j/sin  • | 


2  1  ,,2  S 


.  CO  “  -  ^  V* 
p  4 


HA 


1  2  x** 

j  V  j  +  sin  (t  -  Y) 


This  is  Eq.  (F-3)  in  the  text. 

If  we  further  make  the  assumption  that  (O  »  —  V  ,  then  we  have 

P  4 

5(t)  =  .e/cos  Y  (CO  ^  -  1  +  u  tan  Y)  } - 
•  P  ' 

.  |(exp  (-  vt/2)  sin  Y  cos  co^t 

+  (j  V  sin  Y  -  cos  y]  (sin  ‘*'pt)/‘^pj  +  sin  (t  -  Y) 


At  resonance,  oj  ■*  1, 
P 


tan  Y 


u>p  -1 


+  ®  and  Y  +  7r/2 .  Hence ,  at 


resonance. 


(e/v)  exp  (-  vt/2)  (cos  t  +  i  V  sin  t)  -  cos  t  J 


This  is  Eq.  (F-4)  in  the  text. 


.'•1 


t'lvjure  .  The  right  haiui  siae  of  the  relation 

v'/f  .  1  -  «xr  (  -  4  vt^)(  I  t  4  vtj^] 

plotted  as  a  function  of  the  wave-breaking  time. 
The  dotted  curve  shows  the  wave-breaking  time 
in  the  absence  of  collisions. 

figure  .  The  electron  velocity  at  wave-breakinvj  versus  the 

(scaled)  collision  frequency.  V  -  (::l' )  **  is  tl»e 

o 

electron  velocity  at  wave-break j tig  in  the  case  of 
a  collisionless  plasma. 


'  >»• 


ieriii,ii-i  irieati 


u/feEd/meCoD 
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G.  STUDIES  OF  RANDOM  PERTURBATIONS  ON  PARAMETRIC  EXCITATION  PROCESSES 
One  of  the  more  popular  explanations  for  the  anomalously  strong 
interaction  which  is  observed  between  intense  laser  light  and  plasmas 
involves  the  excitation  of  parametric  instabilities  in  the  plasma.  The 
high  frequency  electromagnetic  field  of  the  incident  laser  beam  couples 
with  the  natural  inodes  of  the  target  plasma  (e.g.,  electron  plasma 
waves  or  ion-acoustic  waves)  cuid  drives  them  unstable.  Such  instabili¬ 
ties  will  eventually  saturate  in  a  turbulent  plasma  state  which  charac¬ 
terizes  the  further  interaction. 

The  two  s inkiest  models  of  such  parametric  process  involve  the 
parametric  excitation  of  simple  harmonic  oscillators: 
i)  Single-mode  parametric  excitation 

^  X(t)  =  0  (G-1) 

dt 

where  the  excitation  enters  through  the  frequency 

^  (1  -  2e  cos  u)  t)  (G-2) 

o  o 

ii)  Coupled-mode  parametric  excitation 

dX 

=  ^1 

,  (G-3) 

d  X 

— 5“  +  2r  —  +  U,  x_(t)  =  X,  z(t)  X,  (t) 

dt2  1  dt  1  2  2  1 

where  Z(t)  =  Z  /2  cos  co  t;  oi,  ,0).  are  natural  oscillation  frequencies 

and  is  the  pump  frequency.  In  this  latter  model,  one  expects 

strong  resonant  energy  transfer  between  the  pump  (e.g.,  incident 

laser  beam)  and  the  natural  oscillations  when  ui  -  oj,  +  u),.  Hence 

0  12 

in  the  usual  analysis,  one  assumes  that 
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a)  a  resonance  condition  is  satisfied:  o)  ~  u, 

o  1  2 

b)  both  the  pump  and  the  damping  are  weak:  e  <<  1,  T  «  co. 

The  general  approach  in  analyzing  parametric  processes  in  plasmas 

is  now  to  cast  the  relevant  equations  describing  the  driven  plasma 
(e.g.  the  hydrodynamics  equations)  into  one  of  these  forms,  and  then 
determine  the  threshold  and  growth  rates  of  ^my  parametric  instekbili- 
ties  which  may  arise. 

In  any  actual  interaction  process,  there  will  be  a  certain  amount 
of  "noise"  which  can  significantly  affect  the  coupling  between  the  pump 
and  the  natural  oscillations.  For  example,  if  the  pump  (laser  beam) 
has  passed  through  a  turbulent  medium,  it  will  be  characterized  by  a 
random  phase  variation.  Furthermore,  most  incident  laser  beams  are 
not  monochromatic,  but  rather  may  exhibit  an  almost  random  frequency 
modulation  as  well  as  a  random  amplitude  variation.  To  investigate  the 
implications  of  such  noise  on  parametric  processes,  we  are  applying 
several  recent  developments  from  the  theory  of  stochastic  differential 
equations. 

The  general  structure  of  such  theories  seeks  to  arrive  at  an 
approximate  equation  for  the  me^ul  solution  <u(t)>  of  a  general 
differential  equation  of  the  form 

3u 

-  L(t)  u(t)  (G-4) 

where  L(t)  is  a  general  stochastic  operator  which  can  be  separated 
into  mean  and  fluctuating  components  as  L(t)  =  <L(t)>  +  6L(t) .  One 
of  the  most  popular  theories  is  based  upon  the  Bourret  approximation 
in  which  <u(t)>  is  described  by  an  approximate  equation  of  the  form 
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t 

dt^  <5L(t)  u(t,t^)  6L(t^)-'«.u(t^)''  (0-5) 

t 

o 

where  the  propagator  u  (t,t  )  is  described  by 

o  o 

9u 

-  ‘■L(t)>  u  (t.t  )  -  0  ,  u(t  ,t  )  -  1  (0-6) 

Jt  o  o  o  o 

By  way  of  example,  if  we  consider  single  mode  parametric  excita¬ 
tion  (ignoring  damping  for  convenience) 


(t)  X(t)  -  0  (0-7) 

dt*" 


and  allow  the  frequency  i*^(t)  to  possess  a  random  component  ^(t) 

i7‘'(t)  -  -  2e  cos  it'^t  c5(t)  (0-6) 


Hence  a  random  frequency'  component  is  seen  to  cause  both  a  damping  as 
wall  as  a  frequency-  shift. 

A  second  example  would  be  coupled-mode  parametric  excitation  with 


a  randomly  modulated  pump 
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Z 

Z(t)  *  C(t)  cos  Cl)  t  (G-10) 

If  we  utilize  the  resonance  conditions  and  the  Bourret  approximation, 
we  find  that  the  random  pump  changes  the  instability  threshold  condi¬ 
tion  from 


Zo^  W 

- 4 - -  (deterministic  pump)  (G-11) 

to 

Zo^  ^1^2 

- »  min  (random  pump)  (G-12) 

that  is,  the  presence  of  noise  raises  the  threshold  by  a  factor 
Aaj/min  (^2' ^2^  where  Au  is  the  spectral  width  of  the  noise.  In  this 
particular  exan^jle,  the  noise  actually  decouples  the  natural  modes. 


31 

H.  DEVELOPMENT  OF  MHD  PLASMA  SIMULATION  MODELS 

In  order  to  understand  the  complex  dynamics  of  nonequilibrium 
plasmas  interacting  with  laser  radiation,  one  must  take  into  accoxint 
the  hydrodynamic  motion  of  these  plasmas  as  they  are  produced,  heated, 
and  expand  into  the  target  chamber.  Of  particular  concern  here  is  the 
theoretical  investigation  of  the  relative  significance  of  various  types 
of  energy  absorption  and  transfer  mechanisms  as  well  as  the  interpreta¬ 
tion  of  experimental  data.  For  that  reason  we  have  placed  a  high 
priority  on  the  development  of  MHD  computer  codes  capable  of  simulating 
exploding  wire  or  pinch  dynamics  when  subjected  to  intense  laser  irradia¬ 
tion.  Most  of  our  effort  in  this  direction  over  the  past  year  has 
involved  the  development  of  nonequilibrium  models  to  describe  the  ioni¬ 
zation  dynamics  and  radiation  emission  from  dense,  metallic  plasmas  of 
low  to  moderate  atomic  number. 

This  work  has  been  stimulated  in  part  by  the  interest  in  explod¬ 
ing  wire  devices  as  intense  X-ray  sources  or  in  using  dense  metallic 
plasmas  to  convert  long  wavelength  laser  radiation  into  X-rays  (using 
line  emission) .  For  elements  of  medium- to- large  atomic  number,  Z,  the 
ion  species  consisting  of  only  one  or  two  inner  shell  electrons  are 
excellent  discrete  X-ray  emitters,  since  their  bound  state  transition 
energies  from  excited  to  ground  level  scales  as  Z  times  the  energies  of 
hydrogen  transitions.  It  is  these  transitions  which  are  the  source  of 
most  of  the  discrete  X-ray  emission  from  exploding  wire  plasmas. 

An  MHD  computer  code  forms  the  basis  for  our  simulation  studies 
in  which  the  macroscopic  plasma  parameters  are  determined  as  a  fvmc- 
tion  of  time  and  radius  in  a  one-dimensional,  cylindrically  symmetric 
configuration  (see  Table  H-l  for  the  relevemt  equations).  The  ionization 
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TABLE  H-1.  MHD  PLASMA  MODEL 
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Here  oi  =  e,  i/  v  is  specific  volume,  u  is  velocity,  q  is  the  usual 
artificial  viscosity  term,  <  is  thermal  conductivity,  tgq  is  equiparti- 
tion  time,  J  is  current  density,  and  £  are  magnetic  and  electric  field 
intensities,  p  is  pressure,  n  is  resistivity,  e  is  internal  energy,  and 
2  is  charge  state.  The  term,  S,  contains  the  plasma  energy  source  and 
loss  terms. 


S 

where  f  is  Joule  heating,  R  is  radiation  emission,  dEj/dt  is  energy 
lost  due  to  ionization  of  the  plasma,  and  Qj  is  the  laser  energy  addi¬ 
tion  term.  The  expressions  for  the  thermal  conductivities  are  taken  from 
Braginskii,  while  the  resistivity  and  equipartition  time  are  taken  from 
Spitzer. 
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&nd  radiation  dyncimxcs  of  dense,  magnetohydrodynamic  plasmas  has  been 
modeled  in  a  fashion  suitable  for  implementation  in  MHD  computer  simu¬ 
lation  codes.  It  has  been  shown  that  more  primitive  LTE  models  such 
as  those  based  on  the  Saha  equation  are  inadequate  to  describe  the 
rapidly  varying  temperature  and  density  regimes  characterizing  many 
such  plasmas.  Both  the  assumptions  of  thermodynamic  equilibrium  in  the 
plasma  and  the  neglect  of  radiative  recombination  in  determining  the 
effective  plasma  charge  state  and  density  distribution  of  ions  cunong 
the  various  states  are  invalid  for  these  plasmas.  By  direct  comparison 
of  theoretical  calculations  based  upon  the  LTE  model  with  those  obtained 
from  the  more  sophisticated  model  based  i^on  time  dependent  rate  equa¬ 
tions  employing  atomic  cross-sectional  data,  and  also  by  comparison  of 
predicted  line  ratios  with  actual  experimental  data,  we  have  shown 
that  the  latter  model  is  necessary  for  an  accurate  estimate  of  plasma 
ionization  dynamics. 

Detailed  collisional-radiativ'e  models  are  developed  which  directly 
solve  the  time-dependent  rate  equations  characterizing  atomic  processes 
along  with  those  equations  characterizing  the  hydrodynamic  motion  of 
the  plasma.  These  models  are  applied  to  analyze  high  density  helium 
2-pinch  and  lithium  exploding  wire  plasmas,  and  they  are  found  to  yield 
results  which  compare  quite  favorably  with  experimental  data. 

A  theoretical  model  of  the  magnetohydrodynamic  (MHD)  motion  and 
radiation  emission  characterizing  dense,  laser-heated  exploding  wire 
plasmas  has  been  developed  which  solves  (numerically)  the  set  of  MHD 
equations  coupled  with  a  set  of  atomic  rate  equations  describing  the 
ionization  emd  radiation  dynamics  of  the  plasma.  The  model  has  been 
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applied  to  study  the  use  of  a  high  powered  laser  to  ionize  and  heat  an 
aluminum  exploding  wire  plasma.  Particular  attention  has  been 
directed  towards  the  radiation  emission  from  such  a  plasma  which 
includes  free-free  (bremsstrahlung) ,  free-bound  (radiative  recombina¬ 
tion)  ,  and  bound-bound  (line  emission)  contributions.  Included  in 
the  line  emission  are  hard  .X-rays  produced  by  excited  level-to-ground 
state  transitions  in  the  helium-li)te  and  hydrogen-li)ce  ion  species. 

Laser  absorption  efficiency  and  laser-to-X-ray  conversion  effi¬ 
ciencies  have  been  studied  as  functions  of  laser  pulse  width,  laser 
intensity,  and  initial  wire  diameter.  It  was  found  that  the  X-ray 
emission  depends  critically  on  the  laser  intensity,  the  plasma  density 
profile  and  hydrodynamic  motion  at  the  time  of  laser  interaction,  and 
the  strength  of  the  plasma  current  and  its  subsequent  self-consistent, 
confining  magnetic  field. 

In  addition,  the  iirportance  of  the  exploding  wire  dynamics  in 
relation  to  the  subsequent  radiation  emission  was  demonstrated  by 
comparing  the  data  on  two  experiments  in  which  the  initial  time  of 
laser  incidence  differed.  .More  specifically,  it  was  shown  that  modify¬ 
ing  the  density  profile,  initially,  for  strong  laser  abosrption  does 
not  assure  what  the  laser  energy  addition  will  be  converted  to  X-ray 
emission  with  maximum  efficiency.  A  careful  balance  of  the  hydrody¬ 
namic  motion  is  required  for  pealc  conversion;  hence,  care  must  be 
taJcen  when  choosing  the  value  of  the  maocimum  current  and  its  frequency 
used  to  explode  the  wire.  Even  then,  the  choice  of  the  appropriate 
pinch  and  the  particular  time  before  that  pinch  at  which  to  trigger  the 
laser  may  critically  determine  the  laser-to-X-ray  conversion  efficiency. 


It  appears  that  the  confining  magnetic  fields  might  be  exploited  to 
contain  the  plasma  at  constant  densities  while  laser  heating  takes 
place.  The  higher  fields  corresponding  to  the  later  time  for  laser 
incidence  (hence,  higher  current  from  the  sine-wave)  in  the  1  mil  wire 
simulations,  no  doubt,  played  a  significant  role  in  the  higher  efficiency 
of  X-ray  production  obtained  in  that  study.  A  modified  current  profile 
that  increases  dramatically  in  a  step  or  ramp  function  after  lase  r  heat¬ 
ing  begins  might  produce  this  confinement. 


I.  EXTERNAL  CIRCUIT  EFFECTS  ON  PLASMA  DYNAMICS 

Our  hydrodynamic  simulation  models  were  modified  to  include  a 
self-consistent  treatment  of  the  coupling  of  the  plasma  to  an  external 
circuit.  Here  we  followed  closely  the  procedure  used  in  the  WHYRAC 
code  developed  at  the  Naval  Research  Laboratory.  This  treatment  models 
the  external  circuit  as  a  driving  voltage  into  tin  affective  lumped  cir¬ 
cuit  characterized  by  resistance  and  inductance.  The  plasma  is  then 
treated  in  an  exact  fashion  as  a  distributed  circuit  element  using  the 
MHD  model.  .An  iterative  procedure  is  employed  to  solve  the  coupled 
Circuit  equations. 

We  have  compared  MHD  calculations  of  the  dyncimics  of  a  fast 
D-pinch  helium  plasma  both  with  and  without  coupling  to  an  external 
circuit  in  Figures  1  and  2.  In  the  latter  calculation  the  total 
current  through  the  plasma  is  specified  as  a  sin  shape  in  time. 

Although  several  detailed  features  of  the  plasma  dynamics  are  modified 
slightly  by  the  inclusion  of  external  circuit  effects  (notably  pinch 
time,  pecik  density,  and  peak  temperature),  the  general  features  remain 


unchanged. 
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TRANSVERSE  EFFECTS  ON  LASER  BEAM  PROPAGATION  IN  PLASMAS 


The  laser-plasma  interaction  experiments  performed  in  our  labora¬ 
tory  indicate  that  an  intense  CO^  laser  beam  can  actually  penetrate 
through  an  over-dense  plasma  column.  A  theory  based  on  nonlinear 
optical  effects  due  to  nonuniform  transverse  temperature  and  density 
profiles  across  the  beam  appeeurs  to  explain  many  of  the  experimentally 
observed  characteristics  of  this  phenomenon.  However  we  also  ^ulalyzed 
this  process  from  a  somewhat  different  perspective  by  making  use  of  a 
plasma  hydrodynamics  computer  code  equipped  with  a  laser  light  propaga¬ 
tion  and  interaction  package. 

More  specifically,  we  modified  the  laser-plasma  hydrodynamics 
computer  code  PHD  to  ^malyze  laser  beam  propagation  through  an  over- 
dense  plasma.  The  PHD  code  solves  the  one  fluid,  two  ten^erature 
equations  of  hydrodynamics  in  one-dimensional  (sleib,  cylindrical,  and 
spherical)  geometry  in  Lagrangian  coordinates.  The  momentum  equation 
is  solved  using  an  explicit  finite  difference  technique  with  the  addi¬ 
tion  of  artificial  viscosity.  The  energy  equations  include  flux- 
limited  thermal  diffusion  through  the  electrons  and  ions,  work  on  the 
electrons  and  ions,  collisional  coupling  between  the  electrons  and 
ions,  and  source  therms.  The  electron  source  term  is  due  to  laser 
light  absorption  calculated  using  inverse  bremsstrahlung  coefficients. 
Laser  beam  propagation  is  analyzed  by  solving  the  relevjuit  photon 
transport  equations. 

The  original  version  of  the  PHD  code  was  not  directly  applicable 
to  the  problem  of  interest  because  it  is  restricted  to  one-dimensional 
symmetry  and  does  not  take  account  of  magnetic  fields.  However  we 
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modified  the  code  to  include  an  approximate  treatment  of  these  effects. 
In  particular,  we  added  terms  to  account  for  magnetic  pressure  effects. 
We  have  modified  transport  coefficients  (most  notably,  thermal  conduc¬ 
tivity)  to  account  for  magnetic  fields.  Cyclotron  radiation  loss  terms 
are  also  included.  To  account  for  transverse  effects,  we  have  included 
artificial  source  sin)c  terms  in  the  hydrodynamics  equations  to  account 
for  trem averse  heat  conduction  and  mass  flow. 

Our  "1.5-D"  hydrodyn^unic  model  explains  several  features  of  the 
experiment.  In  particular,  the  time  required  for  be£un  penetration,  the 
transmitted  laser  pulse  shape,  and  the  basic  characteristics  of  the 
density  profile  are  in  good  agreement  with  the  experimental  observations 
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TRANSPORT  PROCESSES  IN  DENSE  PtASMAS 


It  has  become  apparent  from  both  experimental  measurements  and 
computer  plasma  simulation  studios  that  non-hydrodynamic  transport 
processes  can  play  a  major  role  in  the  dynamics  of  dense  plasmas.  One 
must  provide  a  satisfactory  treatment  of  radiation  transport  effects 
in  any  theoretical  model  of  the  plasma.  Of  particular  concern  hero 
are  X-ray  line  emission,  transport,  and  re-absorption  in  optically 
-  thick  .reqigns  of  ^e  plasma.  Put  situations  can  also  arise,  par¬ 
ticularly  in  laser  heated  plasmas,  when  the  transport  of  non-thermal 
charged  particles  (primarily  electrons)  can  dominate  the  plasma 
behavior. 

Therefore,  we  turned  our  .attention  to  studies  of  how  one  can 
adequately  account  for  such  tr.ansport  processes  within  a  hydrodynamic 
model  of  the  plasma.  Hero,  we  relied  heavily  on  our  experience  with 
very’  similar  problems  which  arise  in  the  study  of  neutron  transport 
in  fission  chain  reaction  systems. 

1 .  Charged  Particle  Tr.insport 

The  accurate  simulation  of  the  transport  of  energetic  (i.e.  super- 
thermaU  cliarqed  particles  is  of  v’ital  importance  to  tlie  success  of  many 
controlled  thermonuclear  fusion  schemes.  For  ex.vnple,  there  lias  been 
great  interest  in  simulating  the  transport  of  supierthermal  electrons 
during  the  implosion  process  in  inertial  confinement  schemes  (notably 
laser  fusion)  as  well  as  investigating  burn-wave  propagation  via  the 
deposition  of  the  1.5  MeV  alpha  particle  energy  in  the  dense  core  of  a 
compressed  microsphere.  In  inertial  confinement  scliemes  using  either 
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If  I 


electrons  or  ions,  there  are  many  questions  to  be  studied  concerning 
the  deposition  of  the  charged  particle  energy  in  the  target,  the  propa¬ 
gation  of  beams  through  space-charge  neutralising  media,  and  the  like, 
that  likewise  require  an  accurate  simulation  of  energetic  charged 
particle  transport.  Charged  particle  transport  also  plays  a  funda¬ 
mental  role  in  areas  of  plasma  physics.  For  example,  there  is  con¬ 
siderable  intere-t  in  obtaining  solutions  of  the  spatially  inhomo¬ 
geneous  Fokker-Planck  equation  to  study  the  energy  deposition  from 
neutral  beam  injection  into  Tokamak  devices  where  the  ionized  beam 
particles  are  either  directly  injected  or  scattered  into  trapped  orbits. 

The  interaction  mechanism  of  most  concern  in  the  transport  of 
energetic  charged  particles  in  plasmas  involves  small  angle  scattering 
collisions . 

The  most  common  description  of  the  cumulative  effect  of  many 
random,  small  angle  collisions  on  a  distribution  of  particles  is  pro¬ 
vided  by  the  Fokker-Planck  equation.  Since  Fokker-Planck  equation  is 
commonly  used  as  the  basis  for  most  numerical  studies  of  charged 
particle  transport  in  plasmas,  a  variety  of  methods  have  been  developed 
to  solve  it.  In  much  of  the  early  work  such  as  that  developed  for 
modeling  magnetic  mirror  devices  the  spatial  dependence  was  suppressed 
emd  only  the  velocity  space  physics  was  retained.  More  recent  methods 
which  specifically  account  for  spatial  effects  include  modified  multi¬ 
group  flux-limited  diffusion  schemes  and  the  integral  particle  track¬ 
ing  techniques.  However,  most  of  these  schemes  involving  the  introduc¬ 
tion  of  substantial  approximations  into  the  Fokker-Planck  equation  in 
order  to  arrive  at  a  set  of  equations  that  are  amenable  to  numerical 
solutions . 
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We  have  developed  a  direct  numerical  approach  which  solves  the 
Fokker-Planck  equation  in  full  generality  including  time,  space,  and 
velocity  dependence.  This  scheme  can  easily  be  generalized  to  include 
more  sophisticated  physics  in  the  future.  To  this  end,  we  relied 
on  the  discrete  ordinates  method  since  it  is  the  most  highly  developed 
numerical  technique  for  obtaining  deterministic  solutions  to  the 
transport  equation. 

The  use  of  the  discrete  ordinate  method  for  solving  the  spatially- 
dependent  Fokker-Planck  equation  is  particularly  attractive  since  it 
has  been  highly  developed  for  the  solution  of  radiation  transport  prob¬ 
lems.  Indeed,  there  are  a  veuriety  of  sophisticated  and  efficient  one 
and  two  dimensional  con^iuter  codes  available  that  have  planar,  cylin¬ 
drical,  and  spherical  geometry  options  (as  well  as  toroidal  and 
triangular  meshes  for  non-orthogonal  geometries)  and  allow  for  a  variety 
of  boundary  and  source  conditions.  Therefore,  the  successful  incorpora¬ 
tion  of  the  Fokker-Planck  collision  term  into  the  multigroup  discrete 
ordinates  formalism  provides  for  an  easy  extension  to  a  variety  of 
geometries  and  source  configurations. 

We  implemented  our  numerical  procedure  for  solving  the  Fokker- 
Planck  equation  by  modifying  a  versatile  time-dependent  discrete 
ordinates  code,  TIMEX,  to  allow  for  the  Fo)cker-Planck  collision  term. 
TIHEX  is  a  time-dependent  particle  transport  code  originally  develops 
for  the  study  of  neutron  transport.  It  utilizes  a  discrete  ordinates 
treatment  of  the  particle  direction  (angle),  a  discontinuous,  linear 
finite  element  treatment  of  the  spatial  variable,  a  multigroup  treat¬ 
ment  of  the  energy  variable,  and  an  explicit  time  differencing  scheme. 
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The  suitability  of  these  features  for  the  analysis  of  charged  particle 
transport  will  become  apparent  as  we  proceed  with  our  development.  We 
refer  to  the  modifications  of  this  code  as  TIMEX-FP.  Typical  results 
cure  provided  in  Fig.  2. 

As  far  as  we  know,  the  TIMEX-FP  code  represents  the  most  con^lete 
attempt  at  solving  the  space,  time,  angle,  and  energy-dependent  Fokker- 
Planck  equation  to  date.  We  have  attempted  to  retain  as  much  of  the 
original  physics  content  of  the  equation  as  possible  and  have  indicated 
that  genearlizations  to  include  nonlinearities  or  spatial  inhomogeneities 
should  be  easily  accomplished  within  the  established  formalism.  More¬ 
over,  the  successful  application  of  established  multigroup  discrete 
ordinates  techniques  has  obviated  the  need  for  independent  and  redundant 
development  of  the  large  variety  of  geometyr,  boundary,  and  source  con¬ 
figurations  that  are  immediately  available  in  a  large  number  of  stock 
production  discrete  ordinates  codes. 

It  has  been  demonstrated  that  the  discrete  ordinates  method  is 
capable  of  accurately  modeling  the  fiill  range  of  the  Fokker-Planck 
collisional  physics  without  recourse  to  any  significant  approximations. 
Therefore,  it  should  provide  a  useful  tool  for  studying  the  details  of 
many  plasma  physics  problems  where,  for  example,  a  knowledge  of  the 
time  dependent  evolution  of  the  single  particle  distribution  function 
or  the  spatial  dependence  of  the  energy  deposition  profile  is  desired. 


Figure  2 
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standard  nmnerical  approaches  to  radiation  transport  such  as  ray 
tracking  and  the  variable  Eddington  method  are  quite  useful  for  time- 
independent  calculations,  but  become  rather  expensive  when  coupled 
to  time- dependent  hydrodyncunic  calculations. 

Therefore  much  of  our  effort  this  past  year  has  been  directed  at 
developing  approximate  radiation  transport  models  more  compatible  with 
implementation  in  a  fast-running  hydrodynamics  code.  Two  approaches 
have  been  studied  thus  far.  The  first  is  based  on  the  frequency  broad¬ 
ening  of  an  emission  line  due  to  Doppler  broadening  during  a  series 
of  eibsorptions  and  re-emissions.  The  second  approach  couples  escape 
and  transmission  probability  techniques  for  optically  thin  regions  and 
diffusion  (Eddington)  theory  for  the  optically  thick  regions. 

The  frequency  broadening  treatment  has  already  been  implemented 
in  the  hydrodynamics  simulation  code.  The  escape  probability /diffusion 
model  has  been  applied  thus  far  only  to  a  few  regions  modeled  problems. 
The  emission  predictions  of  these  two  approaches  will  be  compared  to 
the  results  of  a  variable  Eddington  calculation  to  determine  the 
validity  of  each  model. 

1.  Frequency-Broadening  Model 

The  basic  idea  inherent  in  this  model  is  that  the  line  profile 
widens  to  a  width,  x^,  where  the  optical  depth,  equals  one: 

T,  ”  1  *  exp  (-  X,  *■) 

1  o  1 

where  t  is  the  optical  depth  at  line  center.  The  photons  then  escape 


from  the  wings  of  the  line. 
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The  time  that  it  takes  for  a  photon  to  escape  is  made  up  of  two 
parts:  the  time  for  the  radiation  field  to  equilibrate  (time  for  the 
profile  to  widen  to  and  the  escape  time  (time  for  the  photon  to 
escape  after  the  field  has  equilibrated) .  The  number  of  collisions 
(absorption/re-emission)  that  a  photon  makes  before  escaping  can  also 
be  written  as  the  sum  of  the  number  of  collisions  during  each  time 
interval.  Once  the  number  of  collisions  is  known,  the  probability  that 
a  photon  is  quenched  (photon  energy  removed  from  field  in  a  non- 
radiative  process)  is  given  by 


(1  -  P  ) 

q 


where  P^  is  the  quenching  probability  per  collision.  The  equilibra¬ 
tion  and  escape  times  and  the  number  of  collisions  are  given  by 


t  -  (T  -  1)  t,  ,  t  -  2T  (In  7  )  t. 
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and  t^  is  the  mean  time  between  collisions.  Using  the  quenching 
probabilities,  quenching  rates  are  calculated  for  each  region,  and 
the  resulting  decrease  in  the  radiation  field  is  obtained.  The  net 
energy  loss  due  to  radiation  is  then  used  in  the  electron  temperature 
calculation  in  the  hydrodynamic  simulation. 

2.  Escape  Probability/Dif fusion  Model 

If  a  zone  is  optically  thin,  first  and  second  flight  escape 
probabilities  suffice  to  describe  photon  transport.  The  first  flight 
escape  probabilities  for  the  central  zone  can  be  calculated  in  infinite, 
cylindrical  geometry,  while  for  the  other  regions,  we  can  employ  either 


cylindrical  shell  or  slab  geometries  for  calculating  escape  from  the 
inner  and  outer  surface  areas.  If  all  zones  are  thin,  these  probidoil- 
ities  alone  are  used  to  calculate  total  quenching  and  emission  prob¬ 
abilities  . 

If  a  zone  is  optically  thick,  the  diffusion  equation  with 
appropriate  absorption  and  diffusion  coefficients  is  solved  using 
albedo  boundary  conditions.  The  absorption  rate  within  the  thick  zone 
as  well  as  emission  at  each  surface  is  then  calculated.  Using  these 
quantities  and  the  escape  and  tramsmission  probabilities  of  the  sur¬ 
rounding  thin  zones,  the  total  absorption  and  emission  probabilities 
are  calculated.  These  probabilities  are  presently  being  included  in 


our  MHD  simulation  models. 
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